Controlling complex networks: How much energy is needed? 
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The outstanding problem of controlhng complex networks is relevant to many areas of science and 
engineering, and has the potential to generate technological breakthroughs as well. We address the 
physically important issue of the energy required for achieving control by deriving and validating 
scaling laws for the lower and upper energy bounds. These bounds represent a reasonable estimate 
of the energy cost associated with control, and provide a step forward from the current research on 
controllability toward ultimate control of complex networked dynamical systems. 
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Complex networks are ubiquitous in natural, social, 
and man-made systems, such as gene regulatory net- 
works, social networks, mobile sensor networks and so 
on A network is composed of nodes and edges. The 
nodes represent individual units (e.g., genes, persons, 
sensors) and the edges represent connections or inter- 
actions between the nodes. The state of a node (e.g., 
protein being expressed, opinion of a person, position of 
a sensor) normally evolves over time. And the evolution 
depends not only on the node's intrinsic dynamics but 
also on the couplings with its nearest neighbors (H . 

On one hand, the couplings between nodes increase 
the complexity of collective behaviors, which stimulates 
much interest of modeling, analyzing, and predicting dy- 
namical processes on complex networks Q . On the other 
hand, one may utilize the couplings to control a whole 
network, i.e., steering a network from any initial state 
(vector) to a desired final state, by driving only a few 
suitable nodes with external signals. In this direction 
there are good attempts recently from physics f3-[l0| , bi- 
ology [ill, Il3 ^^"^ engineering [l3l - [l6j research communi- 
ties. Among others, Liu et al. studied the controllability 
of various real- world networks, i.e., the ability to steer a 
complex network as measured by the minimum number 
of driver nodes. A main result was that the number of 
driver nodes required for full control is determined by the 
network's degree distribution Q. Issues such as achieving 
control by using only one controller [l^ and making 
structural perturbations to the network to minimize the 
number of control inputs [loj have also been addressed. 

When control a complex network, an important and 
unavoidable issue is the cost of control. For instance, in 
order to control a social network some efforts has to be 
devoted to change a few individuals' opinions, while to 
control an electronic or a mechanical network, some en- 
ergy has to be consumed to drive a few elements. Even 
if a network is controllable in principle, it may not be 
controllable in practice if it costs an infinite amount of 



energy or if it requires too much time to achieve the con- 
trol. In this Letter, we address this outstanding issue of 
energy cost, i.e., the amount of efforts or energy that are 
necessary to produce external signals for steering a com- 
plex network, and focus on its lower and upper bounds. 
Suppose a complex network is deemed to be controlled 
to a desired state in finite time T/, our main results [see 
Eqs. d?]) and ([5])] show the scaling laws of the energy cost 
bounds with the control time Tf in two different regimes 
separated by the characteristic time. The results give 
faithful estimates for the required energy and thus can 
provide significant insights into bridging network control- 
lability with actual control. 

To be able to analyze the energy cost, we study linear 
networked systems subject to control inputs. This is the 
currently standard framework, upon which the network 
controllability analysis is built 0, Isl-fiol [l6| . A typical 
system of N nodes and M controllers can be written as 



±t = Axi + But, 



(1) 



where xt = [xi{t), X2{t), . . . , XNit)]"^ is the state vector of 
nodes, Uf = [ui{t),U2{t), . . . ,wj\/(t)]'^ is the input vector 
of external signals, B = {bim} is the N x M input ma- 
trix with bim = 1 if controller m connects to node i and 
bim = otherwise, A ~ {o-ij} is the weighted network's 
adjacency matrix including linear nodal dynamics {an}. 

The typical situation of controlling a complex dynami- 
cal network can be characterized as using external signals 
Ut to direct the system Eq. ([1]) from an arbitrary initial 
state Xq toward an arbitrary desired state ktj, in the time 
interval t G [0,TV]. Assuming that the networked system 
is controllable 0, [13] , our goal is to obtain analytic es- 
timate of the energy cost required for achieving control, 
which is defined as [H £{Tf) = J^^ WutW^dt. Gener- 
ally, an infinite number of possibilities exist for choosing 
the control input Uf to steer the system Eq. ([T]) from 
xq to xx, . Of all the possible inputs, the optimal control 



input is given by Ut 



B^eA 



which min- 



2 



imizcs the energy cost [l^, [T^ . The corresponding min- 
imized energy cost is then £{Tf) = v^^'W^^wtj , where 



e^'BRTe^ *dt and vt, = xt, - e 



AT, ^ 



denotes the difference vector between the desired state 
under control and the final state during fi'ee evolution. 
For convenience, we set the origin as the desired state 
= and rewrite the energy cost as 



xo, 



(2) 



where H(T^) = e~ ^Wj'j.e" f is the symmetric 
Gramian matrix [31 . When the system is controllable, 
H is positive-definite (PD), otherwise it is non-invertible. 
In the following wc focus on the normalized energy cost 



EiTf)^£iTf)/\\^,\ 



H 



T 

Xo Xo 



(3) 



When Xo is parallel to the direction of one of H's eigen- 
vectors, the corresponding inverse of the eigenvalue has 
the physical meaning of normalized energy cost asso- 
ciated with controlling the system along the particular 
eigendirection. Using the Rayleigh-Ritz theorem [2^ , we 
can bound the normalized energy cost as 



1 



Vn 



E„un < E{Tf) < En 



1 



(4) 



where Tymax and ?7min are the maximal and minimal eigen- 
values of the PD matrix H, respectively. 

To proceed, we focus on the lower and upper bounds 
of normalized energy cost for the case of single-node con- 
trol. To analytically calculate the quantities 1 /?7max and 
l/??min, for weighted undirected networks, we decompose 
the matrix A in terms of its eigenvectors as A = VSV"^, 
where V is the orthonormal eigenvector matrix that sat- 
isfies VV"^ = V^V = I, S = diag{Ai, A2, . . . , Aat} with 
descending order Ai > A2 > . . . > A^v. We thus have 
e^* = * = Ve^*V'^. Substituting these expressions 
into the Gramian matrix and noting that V is time- 
independent, we have 

H = Ve"^^^' ( / ' e^*V'^BB'^Ve^*dOe~^'^' V"^. (5) 



Denoting the only node under direct control as c, we have 
that B is an X 1 matrix, of which all elements are zeros 
except the cth element, which is one. After some amount 
of algebra, wc obtain 



K,, =Y.Y1 (1 - e-(^°+^^)^^) , (6) 

a=l /3=1 " ^ 

where the Roman letters i,j,c are node indices in the real 
space while the Greek letters a, (3 are running indices in 
the eigenspace. 

To carry the analysis further, we note that there are 
two distinct regimes in terms of the control time Tf. 



In the small Tf regime where Tf <^ l/|Aa -f A^|, we 
can expand e~(-^°+-^'^)'^^' « 1 - (Aq + Xp)Tf and obtain 



TcSicScj- In this 



case, we have Hij « for all i and j except Hcc ~ Tf so 
that the maximal eigenvalue of matrix H can be approx- 
imated as Tf. Consequently, for the small Tc regime, we 
have Emin = l/77max ~ ^/Tf, rcgardlcss of the form of the 
matrix A and of the value of c. In contrast, in the large 
Tf regime characterized by Tf 3> 1/|Aq -f A/3I, we can 
approximate the maximal eigenvalue of H by its trace, 
which has been numerically verified: 77m. 



2A„ 



1-e 



If A is PD, the term e'^^"^/ 



vanishes for large Tf. Wc thus have E„ 



= i/[(A 



A'^)-i]cc- Note that, since the matrix A is independent 
of Tf, the factor 1/[(A + A'^)^^]cc is time-independent 
too. This means that, when A is PD, the lower bound 
of the energy cost converges to a constant value for large 
Tf. If A is not PD, i.e., at least one of A's eigenvalues is 
negative, the most negative eigenvalue Xn will dominate 
the behavior of H: H,, « fl 



'-n ~ 2Aa 

g-2AjvT/^ As a result, the maximal eigenvalue of H 
grows exponentially with Tf. yymax ^ ^-^XnTf tYiat 
Emin ~ (."^^NTf ^ Since Aat < 0, the lower bound of 
the energy cost vanishes exponentially with the control 
time Tf. In the borderline case where A is semi PD, 
i.e., Aq > for a = 1,2,...,A^ — 1 and Xjq = 0, 
the behavior of H can be characterized as: Hj 

-1 



•-■I] 



limA.^o ^-"2^" fl - e-^^-^f) ^ T; 



Our theoretical estimates for the lower bound -Emin 
the energy cost can be summarized as 



of 



En 



1 

[(A+AT)- 



-oxp(2A2vT/) 



small Tf 

large T/, A is PD 
large Tf, A is 



(7) 



PD 



not PD 



Numerical support for Eq. ([7]) is shown in Fig. [TJ We 
use scale-free networks generated by the Barabasi- Albert 
(BA) model [2lj and Erdos-Reyni (ER) type of random 
networks [2^. The link weights are randomly generated 
from the uniform interval [0.5, 1.5]. The linear nodal dy- 
namics arc set as an — —{a + Si) where Si = j^i "^y 
is the strength of node i, and a is such a tunable param- 
eter that one can conveniently change A between pos- 
itive and negative definite. Wc note that other node- 
dependent settings of an will not affect our results. Use 
the method proposed in Q one can find the weighted 
network is controllable, except some pathological link- 
weights sets of measure zero, by any single driver node. 
We numerically compute the lower bound according to 
Eqs. ^ and ©. From Figs. [T^a) and inset of[Hb), we 
see that, for the small Tf regime, E^iin decays as a power 
law regardless of A and c, agreeing with our the- 
oretical result. In the large Tf regime, the behavior of 
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FIG. 1: (color online). Lower bound of the energy cost 
l/rjm!i^ versus the control time Tf. All networks are 
weighted BA scale-free networks except one weighted ER ran- 
dom network in (a), with the same size — 500 and (s) = 20. 
The Sc denotes the strength of directly controlled node. In (a) 
a — —150 which makes A PD. In (b) a = —50 thus A is not 
PD. The dashed line in the semi-log plot in (b) has a slope 
2Ajv. The symbols represent the same quantities calculated 
numerically and the solid lines represent the results from the 
estimation jymax ~ Tr[H]. 



iJmin is determined by the signs of the eigenvalues of A. 
In particular, if the eigenvalues are all positive, the dy- 
namics in the absence of control, i.e., Xt = Axt, will 
force the nodal states to depart away from the zero state. 
Thus, even given sufficiently large time, one has to con- 
sume some amount of energy to steer the nodes back. As 
shown in Fig. [TJa) , E'min converges to a constant value 
as Tf is increased, which agrees with our predicted value 
1/[(A -I- A^)~^]cc- In contrast, if A is not PD, E'^in 
vanishes exponentially, as shown in Fig. [Ijb) . The corre- 
sponding exponent is 2AAr, which is consistent with our 
theoretical estimate in Eq. ([7]) as well. 

We now turn to the upper bound of the energy cost 
Eniax = ^/Vmin- As indicated by Eq. ([5]), most elements 
of the matrix H are small, especially for the small Tr 
regime. Consequently, H is generally ill-conditioned [2^ 
and its minimal eigenvalue is typically very small (though 
positive). Thus, to control a large-size network, E'max can 
be very large. The underlying physical reason is that, 
when only one node is subject to control, the effect on 
other nodes will not be direct but instead will be indirect 
through various paths on the network. The end result is 
that we need to steer the whole system in the state space 
by following highly circuitous, though smooth, routes 
[17| , a process that requires a large amount of energy. 

Typical results computed from Eqs. and ([S]) are 
shown in Figs. [5Ja-c). For small Tf, the upper bound 
Emax exhibits power-law decay, similar to the behavior 
of the lower bound, but the decay exponent for -Emax as- 
sumes a much larger value that is independent of a and 
c [see Fig. [D^d)]. For large Tf, Emax will converge to 
a constant value if A is not negative definite (ND), or 
will vanish exponentially if A is ND. The corresponding 
exponent is given by 2Ai, where Ai is the least nega- 
tive eigenvalue of A, as shown in Fig. [D^e). This is due 
to the fact that, in the large Tf limit, the behavior of 
H~j^ is dominated by the mode with the least negative 
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FIG. 2: (color online). Upper bound of the control energy cost 
Sniax = 1 /rj-mi-n for a Weighted BA network with 20 nodes. In 
(a), a = 2 thus A is ND. In (b), a = -5. In (c), a = -20 
so that A is PD. In (a-c), (•) represent the upper bound 
Smax while (-4) represent the corresponding lower bound E-mi-n 
(included for comparison). In (d) the decaying behavior of 
-Bmax is shown for different Sc and a values. The dash line 
has a slope —36. In (e) the exponential decay of -Emax for 
large Tf is plotted for different values of a. The slopes of 
dashed lines are 2Ai respectively. In (f) the constant values 
of the energy cost in (c) are shown as a function of |a -f- Sc|. 
The slope of the dashed line is 2. 



eigenvalue Ai, which contributes the slowest increase to 
Hij. As a result, we have Emax ^ ^ 

2Ai 



„■ „ V - ~ e'^^^'^f . In the borderline case, i.e., A 

(l-oxp (-2AiT/)) ' ' 

is semi ND, the upper bound decays according to T^^: 

-Einax limAi-yo (i-cKp\-\xiTj)) ~ ^/~^- ^'^^^ ^ behavior 
in both -Bmax and i?min has been numerically verified [l7| . 
The results for the upper bound can be summarized as: 



T7^ (6* > 1) smaU T 



£(A,c) 



■^oxp (2AiT/) 







large T/, A is not ND 
large T,, A is ^^^^ nd 



(8) 



ND 



where e{A,c) denotes a positive value that depends on 
the matrix A and the controlled node c. For the constant 
value of the lower bound as described in Eq. ([7]) , one may 
approximate 1/[(A -I- A"^)~"'^]cc ~ Socc so that Emin is 
proportional to \a + s^l- However, as shown in Fig.[5Jf), 
there appears no proportional relationship between the 
constant value e(A,c) of -Emax and of the controlled 
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FIG. 3: (color online), (a) The ratio eba/eer for different 
network size A'^. In order to eliminate the effects of nodal 
dynamics and strength, we fix the values of an and (s). The 
results include the ratio for optimal driver node (•<) and the 
ratio of averaging over different driver nodes (•). The error 
bars are caused by different generations of network topology 
and link weights, (b) -Bmin {<, left) and -Emax (•, right) 
versus ric, the number of directly controlled nodes. The dot 
pointed by the arrow corresponds to the node with largest 
degree in the network. 



node. This indicates that directly controlling a node with 
larger degree does not generally result in less energy cost. 

Actually, when the system matrix A is PD and 
the control time T/ ^ oo, Eq. [6] reduces to = 

Ea=i EJ=i ^^^IT+xI^'' ^^^""^ solution of AH°° + 
jjcxj^T _ BgT naturally interpreted as dy- 

namical correlation [2^, between nodes i and j with re- 
spect to controlled (driver) node c. So e(A, c) is the in- 
verse of the smallest eigenvalue of the correlation matrix 
H°°. From this point of view, two indications come out 
immediately: Firstly, to find optimal driver node in a net- 
work, one should consider the node viewing from which 
the rest nodes are most dissimilar. The reason is that, 
controlling a central hub node, though may transmit ex- 
ternal signals fast, can induce star-like structure which 
makes the rest nodes more similar to each other. When 
nodes are more structurally similar, they tend to have 
more similar dynamical correlations with other nodes so 
that the corresponding rows in H°° become more similar. 
As a consequence, the smallest eigenvalue of H°° will be 
less. In other words, we have to consume more energy to 
independently steer similar nodes in order to fully con- 
trol the network. Secondly, for randomized networks, 
the more heterogeneous the node-degrees, the higher the 
en ergy cost of control, on the average (see Section III 
of [131 )■ Take randomized BA and ER networks for ex- 
ample, we compare the values of £(A,c), i.e., £ba and 
£er in Fig. 3(a). It shows that the upper bound of en- 
ergy cost for controlling BA networks is much larger than 
that for controlling ER networks. 

We have also studied the energy cost associated with 
the control scheme proposed in a recent work i.e., 
controlling more than one node by a common controller. 



Fig. ISfb) shows the effect of ric, the number of directly 
controlled nodes, on the energy cost, which reveals that 
controlling more nodes will induce smaller value of the 
lower energy bound. This, however, does not hold for 
the upper bound. In fact, adding a node with large de- 
gree into the directly controlled node-set may drastically 
increase the energy cost. This result is consistent, to a 
certain degree, with that found in Ref. Q which shows 
the driver nodes tend to avoid the high-degree nodes. 

It is noteworthy that our results can be easily gen- 
eralized to weighted directed networks. If a network is 
controllable by one driver node, the eigenvalues of the 
corresponding system matrix A are non-degenerate |16| 
though may be not all real. Thus we have A = VSV~ 
where S = diag{Ai, A2, . . . , Aat} with descending order 
of the real part RcAi > ReA2 > ... > RcA^v. Sim- 



ilarly, 
As a 



consequence, 



N 



V.c.(V~^), 



Eq. 



and e 
El is 



St 



(V-i)Te 
replaced by Hij 

1 _e-(^=+^/3)^/V 



Therefore, the scaling laws in Eqs. [7] and H] keep un- 
changed while the decaying exponents are replaced by 
2ReA7v and 2ReAi respectively. Moreover, for large T/ 
and PD A, the constant in Eq. [7]is still proportional to 
2acc by using first-order approximation in [23| . 

In conclusion, we have reduced the complexity of the 
fundamental problem of control cost from the compli- 
cated and intractable Gramian matrix to the simple sys- 
tem matrix which is directly related to the network struc- 
ture. Our results have revealed that energy cost of con- 
trolling complex networks has different scaling behaviors 
with control time in two time scales, separated by the 



characteristic time. 



and 



for the lower and 



2|RoAjv| 2|RgAi| 

the upper bound respectively. In the small-time regime, 
setting a relatively longer time for control always leads to 
less energy cost. While, in the large-time regime, there 
exists the situation where we cannot reduce the energy 
cost even given much more time. Furthermore, our re- 
sults indicate that the lower (upper) bound of energy 
cost is less when controlling a randomized network with 
heterogeneous (homogeneous) node-degrees. These im- 
plications are important when considering the trade-off 
between the energy cost and the control time, which may 
find applications not only for classical [1, Q but also for 
biological [3 [3] and quantum [2J| networks. Al- 
though we have given some heuristics, a method to choose 
an optimal control node-set for minimizing the energy 
cost is lack, which is a promising future work. 
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I. DECAY BEHAVIORS OF E^i^ AND E^^^ IN 
BORDERLINE CASES 

In the main text we argue that, for large control time 
T/, if A is semi positive definite (PD), the lower bound of 
the energy cost i?niin will decay as TJ^ and, if A is semi 
negative definite (ND) the upper bound of the energy 
cost i?max wiU also decay as TJ^ . To provide numerical 
confirmation for these theoretical results, we consider the 
situation of controlling an undirected network of 20 nodes 
which can be controllable by one single driver node. Just 
as in the main text, let nodal dynamics be an = —{a-\-Si) 
where Si is the strength of node i. Setting a = so that 
A is semi ND, we obtain the decay behavior of i?min, 
as depicted in Fig. [SlTa). However, if a = —14.148, A 
becomes semi PD. In this case, we observe the decay of 
Emax as shown in Fig. ISlf b). The dashed lines in both 
figures have the same slope —1. Thus, in these border- 
line situations, E-^aXn and i?max decay as TJ^ for relatively 
large Ty, which is consistent with our theoretical predic- 
tions. 
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FIG. SI: Power-law decays of _Bmin and iJmax for the border- 
line situations. The dashed hnes in (a) and (b) have the same 
slope —1. 



II. OPTIMAL CONTROL ROUTE 



In the main text, we argue that the upper bound of 
the energy cost associated with controlling a complex 
network can be very large, because even the optimal con- 
trol route of steering the whole network from some ini- 
tial state to the origin (zero state) is in general highly 
circuitous, though smooth. Here wc provide numerical 
result of the optimal route for a simple directed network 



6 



0.8 



0.6 



0.4 



0.2 



0.0 



1 ' 1 ' 1 ' 1 ' 1 

Optimal Control Routes - 


...■<£•••••••.. . 








c ^ 


/ - 


1,1,1, 


1 , 1 



-1.0 -0.5 0.0 0.5 1.0 
X,(t) 



FIG. S2: (color online). Optimal control routes to steer a 
simple directed network from the initial state (1.0,0.5)^ to 
the desired state (0,0)^ in the allowed time range [0,1]. 



used in mi: 
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1 



B 



Figure IS2I shows that, when steering the network from 
Xq = (1.0,0.5)'^ to the desired state x^^ = (0,0)"^ in the 
allowed time range [0, 1], the optimal route is indeed not 
direct (dashed line) but smooth and circuitous (•). For 
undirected networks, we obtain similar highly circuitous 
optimal routes (not shown here). 



III. NOTES ON STRUCTURAL EQUIVALENCE 
OF RANDOMIZED NETWORKS 

As we stated in the main text, the constant value 
e(A, c) of energy cost in Eq. (8) is the inverse of smallest 
eigenvalue of the correlation matrix H°° with elements 
= ELi "^"-Ir+^'f-^ Although it is difficult, 
if not impossible, to estimate the smallest eigenvalue of 
H°° , we can obtain some indications from the viewpoint 
of correlations between nodes. When nodes are more 
structurally similar, they tend to have more similar dy- 
namical correlations with other nodes so that correspond- 
ing columns and rows in H°° become more similar. As 
a consequence, the smallest eigenvalue will be less and 
H°° tend to be ill-conditioned so that the upper bound 
of control energy cost increases drastically. 

Because the energy cost of controlling a network de- 
pends not only on the network topology but also on the 
link weights and nodal dynamics, for the convenience 
of comparing the impacts of network structures on the 
control cost, one may set the networks with equal size, 
edge numbers, but different degree distributions, and set 



all nodal dynamics as a function of node-degrees. At 
this situation, when two nodes are structural equivalent, 
the network cannot be controllable by any other driver 
node. Here structural equivalent is defined for any two 
nodes a and b with the same degree as follows: nodes a 
and b are structural equivalent if and only if there exists 
a nontrivial automorphism / of the network |S2| | such 
that /(a) — b. More apparently, take unweighted undi- 
rected networks for example, denote the neighbor-set of 
any node i by S(i), nodes a and b are structural equiva- 
lent if any of the following conditions arc satisfied: 

1. S{a) = 5(6), or 

2. a and b are connected with each other and (<S(a) — 
{b}) = (Sib) - a), or 

3. Vwi e 5(a), 3v2 S S{b) such that {S{vi) — a) = 
{S{v2) — b) or vi and V2 arc connected with each other 
and (iS(ui) — a — V2) = ('5(w2) — b — vi), or 

4. higher-order neighbors of a and b meet condition 3. 
Therefore, for randomized networks, the probability of 

that any two nodes are structural equivalent can pro- 
vide hints for the energy cost of controlling them: the 
larger the probability, the more the energy needed for 
controlling. We here focus on the probabilities of the 
most possible conditions 1 and 2, and ignore the high- 
order probabilities of conditions 3 and 4. 

We will derive the probability of that two nodes are 
structural equivalent for a randomized network with de- 
gree distribution p{k). For a randomized network, two 
nodes i and j are connected with the probability pij = 
j^y^ , where ki and kj are the degrees of nodes i and j re- 
spectively, (k) is the average degree and N is the network 
size. Thus the probability of two nodes a and b with the 
same degree ko have the same neighbors is 



Pcq 



ko 



ko N-2 

N-2Jl\[{k)N) U ^ -{k)N) 



As the degree ki of node i is generated from 
the distribution p{k) independently, the expec- 
tation value of Pcq can be obtained as {pcq} = 
Iki /fe2 ■ ■ • /fc«_2[PcqP(fci)p(fc2) • • .p{kN-2)]dkidk2 ■ . ■ dkN-2 
where the integrations begin at the smallest degree fcmin 
and end at the largest degree fcmax- Therefore, we have 



fco /• / I, \ 2 ^-2 p 

(peq)=^n (m) p(^^)dh n / 



N 
ko 



(k) 
--Z 



p{ki)dki 



where Z = [j^^^) (fco/iV)2(^-2) . For the condition 2 men- 
tioned above, one can obtain similar result. Thus, when 
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the larger the ex- 



fco < N/2, the larger the factor 
pectation value {pcq}- 

From the result one can expect that controlling ran- 
domized Barabasi- Albert networks need much more en- 
ergy than controlling Erdos-Renyi networks, as shown in 
Fig. 3(a) in the main text. Moreover, for randomized 
scale-free networks with degree distribution p{k) oc k~'^ 
(7 > 2), the smaller the value of 7, the more the energy 
needed for controlling. Especially, if ignore the nodal dy- 
namics, when 7 — > 2 the value of pcq ^ 1 for most of 

degrees fcg because is very large at this situation. 
The ending result is that controlling a scale-free network 
with 7 = 2 needs an infinite amount of energy. In other 
words, the network is uncontrollable at 7 = 2 unless most 
of the nodes are set as driver nodes. This indication is 
qualitatively consistent with the result found in |S3l ]. 



IV. REACHABILITY: FROM INITIAL STATE 
xo = TO DESIRED STATE xt^ / 

In the main text, we consider the case of controlling a 
networked system from an arbitrary state xg 7^ to the 
origin xj'^. = 0, which defines controllability. Here, we 
consider the case of steering the system described in Eq. 
(1) from xq = to ^ 0. This situation is referred to 
as reachability [s3 |. 

To analyze the energy cost associated with reachabil- 
ity, we write down the energy expression in the complete 



form as f (T/) = (x^ - xje^ ^f)W^^{xTf 

where Wt^ = J^' e^*BB'^e^^*dt. Since xq = 0, we 
have £{Tf) — x^^W^^^xt^. For undirected networks, 

we factorize A as A = VSV'^, where V is the or- 
thonormal eigenvector matrix with VV"'" = V"^V = I 
and S = diag{Ai, A2, . . . , Ajv} with descending order: 
Ai > A2 > . . . > Atv- Assume that only the c-th node is 
controlled, we have 



Note that the term in the parenthesis is different from 
that of Eq. (6) in the main text. 

Recall the normalized energy cost E{Tf) = £{Tf)/ \\ 
IP, which satisfies 

1/Cmax = -Bmin < E{Tf) < E'max = 1/Cmin, (S2) 

where ^min and ^max are the minimal and maximal eigen- 
values of the matrix Wy^: , respectively. Following a sim- 
ilar analysis in the main text, we approximate ^max by 
the trace of (verified again numerically, as shown 



in Fig. [S3j) and get that 

N 

6nax«^W^,» [=Tr(WT,)] 



1=1 

N N N 



^ ^ 1/ 1/ 

^ A. + a/' 

a=l ^=1 ' 



N 



l)(yy^oy^ 



i=l 



^ ^ T/ T/ 
^ V V 

^ca yea (^2AoTf _ j^-j 



Q = l 



2Aa 



(S3) 

From Eq. (jSSp . we obtain the behavior of the lower 

bound iJmin = 1/Cmax: 



Small Tf regime: l/^^n 



2Ao 



■(e 



2\a,T, 



l)«l/Eo=lKaV^coT/«l/r/.^ 

Large regime: if A is ND, i.e., all of A's eigenval- 
ues are negative, we have 1/^max ~ — 
-1/[(A + AT) 



T\-ll 



2AQ 

If A is not ND, i.e., at least one 



of A's eigenvalues is positive, -Emin will vanish exponen- 
tially as 1/Cmax ~ e~^^i"^^, where Ai is the most positive 
eigenvalue of A. For the borderline situation that A is 
semi-ND, we can obtain the large-time decay behavior by 
setting Ai -> in Eq. (fS3| . which gives i?min ^ 1/^/- 

The behaviors of the lower bound of the energy cost 
associated with the reachability can then be summarized, 
as follows. 



[(A+A'i')- 



-.-cxp (-2A1T/) 



small Tf 

large T/, A is ND 



(S4) 



> large T/, A is 



semi ND 
not ND 



We take unweighted networks for example, and show 
the numerical results on the behaviors of -Bmin versus the 
allowed control time Tf in Fig. IS3I We have also studied 
the relationship between i?min and the degree kc of the 
directly controlled node, as shown in Fig. IS4I We see 
that, if the matrix A is ND, -Emin associated with reach- 
ability will converge to a constant value. In contrast, if A 
is PD, -Eniin associated with controllability will converge 
to a constant value. Furthermore, we recall that, in the 
context of controllability in the main text, driving a node 
with higher degree will induce smaller Emin- However, in 
the context of reachability treated here, driving a node 
with a higher degree will induce larger Emin , as shown in 
Figs. [Ela) andlSHc). 

Following a similar analysis in the main text, we obtain 
the behaviors of the energy upper bound i?max — 1 /Cmin 
associated with reachability, as follows (numerical results 
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FIG. S3: (color online). In the context of reachability where a 
network is controlled from xq = to the desired state xtj 7^ 0, 
the lower bound of the energy cost -Emin ~ 1/^max versus the 
control time Tf . All networks considered are scale-free except 
the one noted as ER Net in (a), which is an Erdos-Reyni ran- 
dom network. The networks have the same size A'^ — 500 and 
the same average degree (k) = 6. In all figures, kc is the de- 
gree of the directly controlled node, and the symbols represent 
the results of our numerical computation and the correspond- 
ing solid lines are the results of our estimation ^max ~ Tr[W]. 
In (a), the parameter is a = 2, which makes the matrix A 
ND. In (b-d), a = -2 so that A is not ND. The dashed 
lines in the log-log plots in (a) and (c) have the same slope 
— 1, which agrees with our theoretical result T^^ for small 
Tf. The dashed line in the semi-log plot in (d) has the slope 
—3.99 ± 0.01, which corresponds to our theoretical estimate 
g-2AiTj,^ as — 2Ai = 2a — —4.0 in this case (see the text for 
details) . 



not shown here). 



T-^iiy » 1) 
^(A,c) 



small Tf 

large Tf, A is not PD 

semi PD 



(S5) 



~oxp {-2XmT[) 



> large T/, A is 



PD 



where <;(A,c) is a positive value depending on A and c, 
and A TV is the eigenvalue of A with the least positive real 
part. 
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FIG. S4: (color online). (a,b) For the context of reachability, 
the lower bound of the energy cost -Emin = 1/Cmax versus 
the degree kc of the directly controlled node. The scale-free 
networks have the same size A*" = 200 and average degree 
(fc) = 6, and the largest degree is km = 54. For comparison, 
(c) and (d) show the corresponding results for the case of 
controllability as discussed in the main text, i.e., steering the 
network from xo 7^ to the desired state xtj = 0. In (a) and 
(b), we set a = 2 so that A is ND and, hence, following the 
analysis in Sec. IIVI we find that -Emin converges to a constant 
value. In (c) and (d), we set a = —80 so that A becomes 
PD and, hence, following the analysis in the main text, -Emin 
converges to a constant value. While (b) and (d) display the 
relationship between the lower bound of the energy cost and 
the degree kc of the controlled single node for different values 
of the allowed control time Tf, (a) and (c) show the behaviors 
for large time Tf = 10, where the node size represents its 
degree, and the node color represents -Emin when controlling 
the corresponding node only. Note that the bottom curves in 
both (b) and (d) show the behavior -Emin oc |a + fcc|, which 
are consistent with our theoretical result Emi-a ~ |1/(A + 
AT)-1| « 2\a + kc\. 



DERIVATION OF OPTIMAL CONTROL ut 



In this paper we have studied the behaviors 
of the energy cost £{Tf) = J^^ ujutdt associ- 
ated with controlling complex networks while choos- 

W7 



mi 
rT, 



g Ut = B'^e^ ''Wrp^^v^j, where vvt/ 

dt, and ~ XTf ~ e^^fxQ. We state 
in the main text without proof that this form of Uj min- 
imizes the energy cost £{Tf). Although the derivation 
of this statement can be found in books on optimal con- 
trol (e.g., [sH|, among others) and is not the subject of 
our present paper, we would like to include it here for 
readers' convenience. 

System: it = Axt -f But, x(0) = xq, x{Tf ) ~ 
Problem: Choose an optimal Ut out of Ut: [0,T/] — > 

to minimize J = J^^ uJutdt. 

The optimal problem can be solved by using Pon- 
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tryagin's Maximum Principle (PMP). Firstly, define the 
Hamiltonian 



Mt = u^ut + A^Axt + But) 



(S6) 



where At is the vector of Lagrange multipliers. Then, ac- 
cording to the PMP, the optimal control signal Ut obeys 



Axt + But, (S7a) 
(S7b) 



A^At 



ut + B^At. (S7c) 



The solution of Eq. (S7b) is At ~ e ^^*c, where c is 
a vector independent of time t. Substituting it into Eq. 
(S7c), we get 



Ut = -B^e A"*c. 



(S8) 



The solution of Eq. (S7a) under the boundary condition 
Xt=o = xo is 



xt = e'^*(xo + f e ^''Busds). 
Jo 



(S9) 



Substituting Eq. (S8) into Eq. (S9) gives Xt = e^*Xo - 
(/o e^^*~''''BB'^e^^^*^'*^(is)e^^^*c. By a change of vari- 



ables t — s = T and thus dr = —ds in the inte- 
gral, we can rewrite the equation as xt = e''^*Xo — 
(J* e^^BBTeA''^dT)e-^^*c. Denoting 



Wt 



(SIO) 



and using the boundary condition Xt=Tf — ^Tf , we have 

= e'^'^^xo - WTfC-^'^'^'c, (Sll) 

where is symmetric as well as PD if the system is 

controllable. It is also referred to as the Gramian matrix 
of controllability. Thus 



where 



„ATf 



(S12) 



(S13) 



Substituting Eq. (S12) into Eqs. (S8) and (S9), we ob- 
tain 

Ut = BTe^"(^/-*)W-;vT,, (S14a) 
Xt = e^'xo + Wte'^^^^J-'^W^VT^, (S14b) 
SiTf) = Joptimal = Wj;Vt,. (S14c) 
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